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We investigate neutron stars in scalar-tensor theories. We examine their secular stability against 
spherically symmetric perturbations by use of a turning point method. For some choices of the 
OS coupling function contained in the theories, the number of the stable equilibrium solutions changes 

os . and the realized equilibrium solution may change discontinuously as the asymptotic value of the 

scalar field or total baryon number is changed continuously. The behavior of the stable equilibrium 
, solutions is explained by fold and cusp catastrophes. Whether or not the cusp catastrophe appears 

D ■ depends on the choice of the coupling function. These types of catastrophes are structurally stable. 

Recently discovered spontaneous scalarization, which is a nonperturbative strong-field phenomenon 
. due to the presence of the gravitational scalar field, is well described in terms of the cusp catastrophe. 
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■^j" . I. INTRODUCTION 

o ■ 

Scalar-tensor theories are among the generalized theories of gravitation. Brans-Dicke theory || is a member 
00 ' of the scalar-tensor theories. Scalar-tensor theories have recently attracted the attention of many researchers. One 
Q\ , of the reasons is that the unified theories that contain gravity as well as other interactions, such as string theory 

naturally predict the existence of scalar fields that relate to gravity. In the hyperextended inflation model [5|, 
q-i scalar-tensor theories of gravity play an essential role. Moreover, projects of laser interferometric gravitational wave 
observations ||J)| will be soon in practical use, so that high-accuracy tests of the scalar-tensor theories may be 
(2jT)! expected pp|-[l4l- 

Scalar-tensor theories are viable theories of gravity for some choices of the coupling function which is contained in 
. ' the theories. Predictions of these theories in a strong field may be drastically different from those of general relativity. 
Recently, Damour and Esposito-Farese Jl5|,[l6| discovered one example of such phenomena. They showed that, for 
some choices of the coupling function, the configuration of a massive neutron star deviates significantly from that in 
general relativity, even if the post-Newtonian limit of the theory is extremely close to or even agrees with that of 
general relativity. This deviation in a strong field may be easily tested from binary- pulsar timing observations, if it 
exists, because of the extra energy loss by scalar gravitational radiation [jl6). The deviation from general relativity 
can be no longer dealt with as a perturbative effect from general relativity. Damour and Esposito-Farese referred to 
this nonperturbative strong-field effect as "spontaneous scalarization" in analogy to the spontaneous magnetization 
of the ferromagnets. 

In this paper, we investigate spontaneous scalarization in detail with the technique of catastrophe theory. A 
many-parameter version of the turning point method fl^-ffiifl is used as a tool of a stability analysis of equilibrium 
solutions. The stability analysis of boson stars in scalar-tensor gravity via catastrophe theory in the case of one- 
dimensional control space was done in |2l] , ^| . Catastrophe types of neutron star equilibrium solutions are classified 
as fold and cusp catastrophes. The occurrence of the cusp catastrophe depends on the choice of the coupling function. 
Spontaneous scalarization is classified as the cusp catastrophe. From this catastrophic feature, we conclude that 
the stable configuration of the neutron star may change discontinuously as the baryon number of the star or the 
asymptotic value of the scalar field changes continuously. The behavior of the scalar charge around the cusp point is 
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explained by catastrophe theory. For the coupling function considered here, when the asymptotic value of the scalar 
field is such that the theory agrees with general relativity in the post-Newtonian limit, we find the sequence of the 
equilibrium solutions bifurcates to three branches at some critical central density. One branch consists of solutions 
that are identical to neutron stars in general relativity, and the other two consist of solutions that deviate significantly 
from neutron stars in general relativity. The general relativistic branch is secularly unstable in agreement with the 
result obtained by a perturbation study [P3[ , while the non-general-relativistic branches are secularly stable. 

This paper is organized as follows. In Sec. II, we summarize the field equations of scalar-tensor theory and 
the equations determining equilibrium solutions of neutron stars in this gravitational theory. In Sec. Ill, we present 
stability criteria on the grounds of the turning point method. In Sec. IV, the stability criteria are applied to equilibrium 
solutions of neutron stars in scalar-tensor theory and some consequences of catastrophe theory are discussed. Section 
V is devoted to conclusions. We use units in which c = 1. The Greek indices run from to 3. We follow the 
Misner-Thorne- Wheeler p3 sign conventions for curvature quantities. 



II. BASIC EQUATIONS 

Here we consider a class of scalar-tensor theories in which gravity is mediated by not only a metric tensor but also 
a massless scalar field. The action is given by [|| 

I= 16^G~ J ^^( R *~ 2 9*" ( P^ l P^ dix + I ^ m ,A 2 ^)g^}, (2.1) 

where g*^ is the "Einstein" frame metric tensor, <P m denotes matter fields collectively, and G* is some dimensionful 
constant. In this Einstein frame, the Einstein-Hubert term is isolated from other sectors. The "Brans-Dicke" frame 
metric tensor g^ v is related to the Einstein frame metric tensor by the following conformal transformation: 

9nv = A 2 (ip)g* lllJ . (2.2) 

Beca use of the "universal coupling," which is the way of the coupling of the scalar field in the matter sector seen in 
Eq. (2T), a test particle moves on the geodesic of the Brans-Dicke frame metric g^. For this reason the Brans-Dicke 
frame is often called a "physical" frame. The tilde denotes the physical frame quantity. 
In the Einstein frame, the field equations are given by 

G* M „ = 87rG*T*^ + 2 \<p^ip, v - \g*^g'i l3 ^, a ^,fj J , (2.3) 



□*y> = -4vrG»a(^)T», (2.4) 

while the equations of motion for matter are 

V*„T^ = a(ip)T^V^^ip, (2.5) 

where the energy-momentum tensor of the matter, T* v , is defined and related to the physical energy- momentum 
tensor as 

rp^v _ 2 dl m [^ m , A 2 (ip)g* fll/ ] _ A 6 f(p\'T fa ' (2 6) 
\J g* &g*fj,v 

G*^ and □„ are the Einstein tensor and d'Alembertian of g*^, respectively. T* and a(tp) are defined as 

n = v = T r<?*^, (2.7) 

, . dhiA(ip) . . 

, — ■ (2-8) 

dip 



The parameters in the parametrized post-Newtonian framework are given by 

Jag 
l + < 

2(1+^)2 

C = cti = ct 2 = a 3 = 0, (2-11) 



1-7^=7^2- (2-9) 
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where 0Edd and jEdd are the so-called Eddington parameters. We have defined 



a = a(ifi ), 



(2.12) 
(2.13) 



and ipo is the value of the scalar field tp in the spatial asymptotic region. We assume that the cosmological evolution 
of the scalar field is sufficiently slow in comparison with the characteristic time scale of the local gravitational process 
of the isolated object considered here. From this assumption ip is regarded as the cosmological value of the scalar 
field. On the other hand, we can identify the asymptotic value ipo to the value of the scalar field in the matching 
region in the matching approach to the TV-compact-body problem (see Appendix A of ||]). Then, the solar-system 
experimental constraints are [E5| 



and 



lEdd = 0.9996 ±0.0017 



^Edd - lEdd - 3 = -0.0007 ± 0.0010. 



(2.14) 



(2.15) 



We summarize equations for the structure of a relativistic star in scalar-tensor theory, following |15|. We restrict 
ourselves to the static and spherically symmetric case. The metric is given in the following form: 



dsi 



^dt 2 



1 



Mr) 



dr 2 + r 2 {d0 2 +i 



The matter is described as a perfect fluid: i.e., 

Tpv = (p + p)u ll u u + pg^. 
Then, the following equations are obtained: 

p! = AnGyA i p + ^r(r - 2p)^ 2 , 



8irG* 



r - 2p 



7"0 



2// 



r(r — 2p) 



ip' = ip, 

ip' = AttG* - — — [a(p — 3p) + r(p — p)tp] — — — J^ip, 



P = -(p+P) 
P = P(p), 



2p l yr yr r{r-2p) 



where the prime denotes a derivative with respect to r. We use the polytropic equations of state: 

r 



p = nm b + 



Knoirib ( n 



p = Knoirib 



n 



n 

1.66 x 10~ 24 



r - 1 \n Q 
r 



0.1 fm 



-3 



N = 



AirnA^r 2 



(2.16) 
(2.17) 

(2.18) 

(2.19) 
(2.20) 
(2.21) 

(2.22) 
(2.23) 

(2.24) 

(2.25) 

(2.26) 
(2.27) 



where h is the baryon number density in the Brans-Dicke frame. We then take the parameters T = 2.34 and K = 0.0195 
(EOS II of JTH]). Note that the total baryon number is given by 



(2.28) 
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Here we present the method of solving the above equations and obtaining the structure of a neutron star. First the 
initial values of the above set of ordinary differential equations are fixed as 



H(0) = 0, 1/(0) = 0, p(0) = (p c , V(0) = 0, p(0) = po, (2.29) 



and Eqs. (2.1S)-(2.23) are integrated numerically up to the stellar surface at which p = 0. Thereafter the solution 
is matched with the static and spherically symmetric "vacuum" solution, where the term "vacuum" means only the 
absence of matter, i.e., T^ u = 0. This solution is given in Qj. The solutions are parametrized by three parameters b, 
d, and ipo. From the matching conditions at the surface, we can obtain v(r) including the constant term. In order to 
set A(ipo) to unity, we rescale the raw quantities to the renormalized ones as follows: 

Tren — AqT : Pren — AqP, ^ren — ^ Pren = <Pi ^Pren = 

tfioren = fo, a ren = A a, b ren = A b, d ren = A Q d, N ren = A^N. (2.30) 

Then, from the asymptotic properties at spatial infinity for a static and isolated system 

2G„m r „f(%\ , 00 ^ 
9*nv = t)iiv + — <V + ° ) > (2-31) 

r ren \ r ren J 

<P = <Po + ^+o(^-), (2.32) 

' ren \ 1 ren J 

where r]^ is the Minkowskian metric. We call m the Arnowitt-Deser-Misner (ADM) energy and uj the scalar charge [|] . 
b ren and d ren are related to m and u> as 

G*m=^, (2.33) 
G*uj — —d ren . (2.34) 

Hereafter the subscripts "ren" are omitted for simplicity. We use units in which G* = 1. 



III. STABILITY CRITERIA 



In scalar-tensor theory, control parameters of the static, spherically symmetric, and isolated neutron star are not 
only the baryon number N but also the "external field", that is, the asymptotic value of the scalar field ip®. For static 
systems, the partial derivative of m in terms of (po with N constant is given by 



dm 
d<p 



N 



The energy injection of the system by increasing baryons is described as 

-1/2 



j 4nu SpA 3 r 2 (l - Q\ dr = J Aire"' 2 pSnA 4 r 2 (l - ^) ' ' dr, 
where p = dp/dn is the chemical potential. The first law of thermodynamics in an adiabatic process is 



-pd 



(3.1) 



(3.2) 



(3.3) 



From Eqs. ( 2.1£ ), ( 2.2C ), ( |2.22|) , and (3.3), we find that the quantity Ae v l 2 p is constant all over the star. Therefore 
this quan tity can be estimated by its value at the stellar surface. Using this fact, the expression of the energy injection, 
Eq. (3.2), is rewritten as 



A s e u ° /2 ~p s \ 4irA 3 r 2 ( 1 



2/< 



-1/2 



Shdr = A s e"° /2 p s 6N, 



(3.4) 



where the suffix "s" indicates that the quantity is evaluated at the stellar surface r — r s . Therefore, the effective 
chemical potential, p e ff, is given by 
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Me// 



\dN J 



(3.5) 



From the above discussions the variation of to for static systems in a quasistatic process is written in the following 
form: 



5m = —Ludipo + ii e ffSN, 



(3.6) 



where by "quasistatic process" we mean successive changes among the infinitesimally nearby equilibrium solutions. 

Suppose that the isolated neutron star is perturbed slightly with ipo and N constant for some reason other than 
incident waves, while spherical symmetry is preserved. Then the outgoing waves can carry out some positive energy 
to infinity and the system cannot keep its original state if there exists an energetically favorable configuration which is 
infinitesimally deformed from the system with the same tpo and N. Therefore, an equilibrium solution X is secularly 
stable against spherically symmetric (infinitesimal) perturbations if and only if there is no momentarily static and 
spherically symmetric configuration Y which is arbitrarily close to X with the same ipo and N but strictly smaller to. 

In order to examine the stability of the equilibrium solution, we follow the turning point method Ej-pOt|. In the 
present problem, to is a potential function, since the equilibrium solution is a stationary point of to (5m = 0) and the 
stable equilibrium solution is a minimal point of m (5m — and 5 2 m > 0). The asymptotic value of the scalar field 
ipo and baryon number N form a two-dimensional control space. The equilibrium solutions are uniquely parametrized 
by two parameters, i.e., the central value of the scalar field, tp c , and the central total baryonic density, p c . 

We adopt the following stability criteria [ po| : 

(i) The stability of X(ip c ,p c ) can change typically only at a "turning point." Here the "turning point" (y^iPc) i s a 
point where there exists a nontrivial vector (5<p c ,5p c ) such that 



r»r fdN\ . fON\ n 



From Eqs. (3.7) and (3.S), 



d((fc,Pc) 







(3.7) 
(3.8) 

(3.9) 



at the turning point. Therefore the change of stability can be detected as envelopes of a family of curves p c = const 
in the (<po,N) plane. Of course, this is also true for a family of curves <p c = const. 

(ii) In order to specify an unstable branch at the turning point, we draw the sequence of equilibrium solutions in 
the ((po,Lu) plane, maintaining TV constant. Then, as one proceeds along the curve in a counterclockwise direction, a 
branch beyond the turning point is unstable. This is also the case with the curve (N, —p, e ft) with ip Q constant. This 
is a direct consequence of theorem I of j2(J] . 

Here we describe the meaning of criterion (i) in the context of catastrophe theory. We regard the ADM energy to 
as a function of three variables tpo, N, and u>. We take u> as a state variable. We consider an equilibrium space 



M„ 



(<po,N,u) 



dm 



= 



(3.10) 



<Po,N 



and a control space 



R 2 = {(^ ,A0}. 



(3.H) 



We define a catastrophe map 



Xm ■ M m 

(ip ,N,uj) 



R 2 , 



(3.12) 



A point P £ M m is called a singular point of Xm if the Jacobian of Xm vanishes at P. A point Q £ R 2 is called a 
singular value if there is at least one singular point in XmiQ)- A bifurcation set B m C R 2 is a set of singular values. 
At a singular point P £ M m , a vector normal the tangent space of M m , which is 
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d 2 m \ / d 2 m \ 
d(p du) N ' \dNdu;) l 



duj 2 



is parallel to the (po-^ plane. Therefore, the set of singular points, S m c M m , satisfies 

' d 2 m s 



and the bifurcation set S m C R 2 satisfies 

B m = J (po.JV) 



/ <9ra\ 



eta 2 







dm 



(d 2 



<Po,N 



\duj 2 



= 



(3.13) 



(3.14) 



(3.15) 



V>o,N 



The envelopes of the family of the curves p c = const in the (<po, N) plane form a bifurcation set B m of the catastrophe 
map Xm because the Jacobian of x-m vanishes at points on the envelopes. Criterion (i) says that a sequence of the 
equilibrium solutions can change its stability only at the points of the bifurcation set. 

From criteria (i) and (ii), we examine the stability of the equilibrium solutions of neutron stars in scalar-tensor 
theory. From the turning point method alone, however, we cannot say that an equilibrium solution is stable. Therefore 
the stability of an equilibrium solution must be investigated by perturbation study once for all. For this purpose, we 
examine the case in which a(<pa) = and m/r s is sufficiently small (p c is sufficiently small). In this case, there is an 
equilibrium solution that is identical to that in general relativity. For this solution, the second-order variation of m 
by regular, adiabatic, time-symmetric, and spherically symmetric perturbations with (po and N constant is 



S 2 m = general relativistic terms + — 



dre -W2 I 1 _ 2 ± 



-1/2 



d 2 



c, 



where 





rScp, 


dr* = 








V(r) = 


r \ r J 



\ 1/2 

dr. 



M r - n 
r(r — 2/i) 



e v -4 7 r/3o(- / 5 + 3p)e", 



(3.16) 

(3.17) 
(3.18) 

(3.19) 



and see Appendix B of (£7J for the general relativistic terms. The general relativistic part is positive definite if T > T c , 
and r c — > 4/3 in the Newtonian limit |2q] . The second term is positive definite because the eigenvalues of the operator, 
—d 2 /dr 2 + V, are all positive for an arbitrary coupling function A(ip) if m/r s is sufficiently small, which has been 
shown in pS}] . Therefore, the general relativistic equilibrium solution in which the central density is sufficiently small 
is stable for the case a(<po) = if Y > T c ~ 4/3. 



IV. RESULTS 



Hereafter we restrict our attention to the coupling function of the quadratic form 

1 



A(tp) = exp ( -/V 



(4.1) 



For th is mo del, t he solar-system experiments constrain the present cosmological value of the scalar field through Eqs. 
(|2.14[) and fl2.15|) as 



\<Po\ Z 0.032| 



and 



,< , 0.012(1 + /3)- 1 /2| 
~ 1 0.029|1 + /3|- 1 /2| / 



- 1 for $ > -1, 
1 for < -1, 



(4.2) 



(4.3) 



respectively. In particular, if (po = 0, the post-Newtonian limit of this theory agrees completely with that of general 
relativity because a(ipo) = 0. 



G 



A. (3 > -4.35 case 

We present here the results of the case f3 = —4, but the features are basically common to the case /3 > —4.35. 
Figure 1 shows p c = const curves in the (ipo,N) plane, where the equilibrium solutions have been determined in the 
manner described in Sec. II. At a point on an envelope of the family of the curves seen in Fig. 1, the stability of the 
sequence of equilibrium solutions changes. Figure 2 shows the curves (po,u>) with N constant. In Fig. 2, the solid 
lines denote stable branches while the dotted lines denote unstable branches, where stability criteria (i) and (ii) are 
applied. Therefore, in region (A) in Fig. 1, only one stable equilibrium solution exists. For ipo ~ 0, this stable solution 
is identical to that of general relativity. In region (B) in Fig. 1, however, no stable solution exists. This is classified as 
the fold catastrophe in which the control space is two-dimensional. This catastrophe is elementary and structurally 
stable. Hence it is expected that this catastrophe structure is not changed by adding small higher-order terms to the 



exponent of the coupling function (4.1). The potential function m is written locally around point p(ipo p , N p ) on the 
envelope (see Fig. 1) as, for po p > 0, 

m = -(u - uj p f + [B(tp 0p - ipo) + B'{N - N p )](u - lu p ) + m p , (4.4) 

where A, B, and B' are some positive constants. Then the terms in the square brackets cancel out on the envelope, 
and are negative in region (A) and positive in region (B). For po p < 0, replace lo — uj p and po p — po with oj p — u> and 
i po — ipo p , respectively. For simplicity, we describe the behavior of the scalar charge for the case po p > 0. From Eq. 



(4.4), near point p, the scalar charge is given by the roots of the following quadratic equation: 

=A(u J -u J p) 2 + [B(ipop~(po)+B , (N-N p )} = 0. (4.5) 

dUJ Ko,N 

The scalar charge is then given near point p in region (A) by 

uj = ujp ± A~ 1 /' 2 [B((p - pop) + B'(N p - TV)] 1 / 2 , (4.6) 
where the upper sign denotes the stable branch and the lower denotes the unstable one. If a quadratic term in (uj — lu p ) 



was invol ved in Eq. (4.4), the number of the roots of the equation dm/du; — did not change at point p. That is 



why Eq. (4.4) does not contain the quadratic term. The "scalar susceptibility" \v ^ gi ven near point p by 



(^) =±U-^ 2 B[B(p a -p 0p ) + B / (N p -N)}- 1 / 2 . (4.7) 



The bifurcation set B m C R 2 , which is the envelope seen in Fig. 1, is given by 

— = A(« - up) 2 + [B(p 0p - p ) + B'(N - N p )] = 0, (4.8) 



dm 
dui 2 



= 2A{u-u) p )=Q, (4.9) 



B' 

<Po = -^(N - N p ) + p 0p , (4.10) 

near point p. This fold catastrophe appears also in general relativity in which A(p) = 1 identically. In general 
relativity, because of the absence of a gravitational scalar field, the control space is one-dimensional. For j3 = —4, the 
maximum ADM energy is greater than the general relativistic one for cpo ^ 0. This is because, due to the presence 
of the scalar field, the effective gravitational constant becomes smaller and thereby gravity becomes weaker than in 
general relativity. 

B. /3 < -4.35 case 

This case is more interesting than the above case. We present the results of the case (3 — —6. Figure 3 shows 
p c = const curves in the (tpo,N) plane. This figure is very different from Fig. 1. On the envelope of the family of 
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curves, e'dcbab'c'de, the sequence of equilibrium solutions changes its stability. Although there are other envelopes in 
region (B), they have nothing to do with the change of the number of stable equilibrium solutions. Figure 4 shows 
curves (ipo,uj) with N constant. From criteria (i) and (ii), the number of stable equilibrium solutions is as follows: 
In region (A), only one stable equilibrium solution exists. In region (B), two distinct stable equilibrium solutions 
exist. Surprisingly, these stable equilibrium solutions are different even for ipo = from their counterparts in general 
relativity. For ipo — 0, the unstable solution agrees with the stable solution in general relativity. One of the two 
stable equilibrium solutions disappears on the envelope dcbab'c'd. In region (C), no stable equilibrium solution exists. 
Point a is a bifurcation point for ipo = 0. This is seen in Fig. 5 which displays the curves (p c ,m) and (p c ,N) for 
ipo = 0, where the solid lines denote stable branches and the dotted lines denote unstable branches. In this figure, 
two stable branches are degenerate because, for ipo = 0, two stable equilibrium solutions are identical except for the 
sign of the scalar field. The equilibrium solution of the bifurcated stable branches is more compact for smaller mass 
but less compact for larger mass than the general relativistic sequence. Figure 6 shows the equilibrium space M m 
near point a. This type of the catastrophe at point a is classified as the cusp catastrophe in which the control space 
is two-dimensional. The map Xm is a cusp catastrophe map. This catastrophe is elementary and structurally stable, 
which su gges ts that this structure is stable against adding small higher order terms to the exponent of the coupling 
function (4.1). Point a is called a cusp point. 

We restrict our attention to cusp point a(0,N a ). (nibN a ~ 1.24M© for (3 = —6.) The potential function m is 
written around cusp point a as 



C 



D(N — N a ) 2 



U3 



ip uj - 



(4.11) 



where C and D are some positive constants. The reason why the coefficient of ipoiv is determined is that Eq. (3.1) 
holds. This form of expansion agrees with the usual Landau ansatz for a second-order phase transition, which has 
been used to explain spontaneous scalarization by Damour and Esposito-Farese [ pjSf . The scalar charge is given by 
the roots of the following cubic equation: 



f)fn 

) = Cuj" - D(N - N a )u - = 0. 



<Po,N 



From Eq. ( 4.12Q , near cusp point a, the scalar charge is given by 

to = 0, 



(4.12) 



(4.13) 



for N < N a with ip = 0. This is a stable branch. For N > N a with <p = 0, 



±(§) 1/2 (iv~AgV2 





for the stable branches, 
for the unstable branch. 



(4.14) 



At point a the stable equilibrium solution changes continuously, but its derivative with respect to N is discontinuous. 
If Eq. (4.11) involved a cubic term in u>, the number of roots of the equation dm/duj — did change at point a. 
But this catastrophe was classified as the fold type and therefore not the case for point a because of the shape of the 
bifurcation set seen in Fig. 3. That is why Eq. (4.11) does not contain the cubic term. We also note that, for the 



case of two-dimensional control space, the structurally stable catastrophe is classified as either the fold or cusp type 
by Thorn's theorem. Therefore, at point a, a second-order phase transition occurs. If we fix N to N a , 



(4.15) 



This is stable. From Eqs. ( |4.12| )-( |4.14| ), with ipo = near point a, it is derived that the scalar susceptibility \v is 
given by 



hD-^N-Na)- 



for N < N a , 
for N > N a . 



(4.16) 



Near point a in region (A) in Fig. 3, only one real root of the cubic equation (4.12) corresponds to the stable equilibrium 
solution, while, in region (B), the smallest and largest roots of three real roots correspond to the stable equilibrium 
solutions and the intermediate root corresponds to the unstable one. If (fo > 0, the largest root corresponds to the 
globally stable one. If ipo < 0, the smallest root corresponds to the globally stable one. If ipo = 0, the two stable 
equilibrium solutions have identical ADM energies. The bifurcation set B m , which is the envelope b'ab, is given by 
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dm 
dui 

d 2 m 
IhS 2 



= Clu 3 - D(N - N a )uj -<p = 0, (4.17) 

<Po,N 

= 3Clu 2 - D(N - N a ) = 0, (4.18) 

Vo-N 



i.e. 



<P° = ±(^ d ) (^-^) 3/2 , (4.19) 

near point a. The cusp catastrophe has been named after this shape. On the envelope dcbab'c'd except for points a 
and d, one of the two distinct stable equilibrium solutions, the locally but not globally stable one, disappears, and 
hence a first-order phase transition occurs if the system obeys a perfect delay convention. On envelope ede', the stable 
equilibrium solution disappears. The catastrophic feature on the envelopes except for point a is the fold catastrophe 
described in the last subsection. Point d is not a cusp point but the intersection of two folds. 

We should comment that, for the near critical case —4.9 < j3 < —4.35, the behavior of the stable equilibrium 
solutions around the point of the maximum baryon number is somewhat complicated, although the structure of the 
cusp catastrophe at cusp point a is not changed. Figures 7 and 8 show the curves (p c , N) with <po = 0, for (3 — —4.5 
and —4.85, respectively. For —4.8 < (3 < —4.35, the number of stable equilibrium solutions changes as 1, 2 (degenerate 
in Fig. 7), 3, 1, as the control parameter N is increased continuously, as is seen in Fig. 7. For —4.9 < (3 < —4.8, the 
number of stable equilibrium solutions changes as 1, 2 (degenerate in Fig. 8), 3, 2, as N is increased continuously, 
as is seen in Fig. 8. For (3 < —4.6 the maximum ADM energy with po = is greater than that in general relativity, 
while, for (3 > —4.6, it is the same as that in general relativity. 

The Kepler mass, which governs the Newtonian orbital motion of a test body, is not the ADM energy m in general, 
but§ 

1 + a a a A , , 

V = -i ~ 2 TO ' ( 4 ' 2 °) 

where 

d In m u> , , „„ . 

a A = ^—= . (4.21) 

<j(p m 

When we consider the case of ao = /3<^o = 0, the Kepler mass is identical to the ADM energy. Therefore the argument 
above for ipo = is also valid for the Kepler mass. 

Here we present the physical interpretation as to why spontaneous scalarization occurs. In spite of the absence of 
the potential in the Lagrangian, the scalar field <p obtains an effective potential term W(ip) which satisfies 

dW 

— = -4ira(<p)T* (4.22) 



because of the coupling with matter. Note that T* depends on ip. Then, if we consider A(ip) of the form (|4.l| ), 



dV 

— = -47r/fy>T„ (4.23) 
dip 

and if T* = A A {—p + 'dp) is negative, ip = is an unstable stationary point of the effective potential, if (3 < 0. On the 



other hand, the term from the spatial derivative in Eq. (2.4) has a contribution to stabilize the solution. By these 



two competing effects, the stability of the trivial configuration <p = against spontaneous scalarization is governed. 
For a detailed analysis of the stability of the trivial configuration, see p3[ . 

If spontaneous scalarization occurs, the effective gravitational constant, which is A 2 (ip) = exp(/3p 2 ) in the sense of 
the inverse of the Brans-Dicke scalar field, becomes considerably smaller than unity. Thereby the gravitation becomes 
weaker and a considerably larger mass than in general relativity can be supported by the lower matter pressure than 
in general relativity. 
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V. SUMMARY AND DISCUSSIONS 



The behavior of the equilibrium solutions of neutron stars in scalar-tensor theories of gravitation shows a catas- 
trophic feature, which is characterized by a discontinuous change of the system. When we consider a function A(ip) 
of the form A(ip) — exp(^(3(p 2 ), the catastrophe types are classified as fold and cusp catastrophes. The appearance 
of the cusp catastrophe depends on whether (3 > —4.35 or /3 < —4.35. From the fact that those types of catastrophes 
are structurally stable, it is expected that they would be seen in a wide class of coupling functions. 

For (3 > —4.35, the fold catastrophe on the two-dimensional control space does occur. The critical baryon number 
and critical ADM energy depend on tpQ. For a baryon number smaller than the critical one, one stable equilibrium 
solution exists, while, for a baryon number larger than the critical one, no stable equilibrium solution exists. In 
particular, for (po — 0, the stable equilibrium solution is completely identical to that in general relativity. The 
behavior of the scalar charge and scalar susceptibility near the critical baryon number is explained by the form of the 
potential function of the fold catastrophe. 

For j3 < —4.35, the cusp catastrophe does occur while the fold catastrophe also occurs. For (3 < —4.9, there is 
some critical value of the scalar field, <fiQ rlt > 0. If l^oj > i^o"*' there is only one critical number N cntl that depends 
on (pg. For N < J\[ mtl , one stable equilibrium solution exists, while, for N > N cntl , no stable equilibrium solution 
exists. If < \<po\ < v?o"*, there are three critical baryon numbers N critl > N crit2 > N crit3 . For TV < jV crit3 or 
jycnt2 < TV < N cntl , only one stable equilibrium solution exists. For N crlt3 < N < jV c "* 2 , two distinct stable 
equilibrium solutions exist and they do not agree with those in general relativity even for the limit (po — > 0. The 
almost general relativistic branch is unstable for N > N crlt2 . For N > N cntl , however, no stable equilibrium solution 
exists. If ipo — 0, the sequence of equilibrium solutions of neutron stars bifurcates at a point. Beyond this point, the 
general relativistic branch becomes unstable and another two (degenerate) sequences of equilibrium solutions far from 
the general relativistic one are stable. This bifurcation point is a cusp point, and the behavior of the scalar charge and 
scalar susceptibility near the cusp point is explained by the form of the potential function of the cusp catastrophe. At a 
point on the envelopes other than the cusp point, the fold catastrophe occurs. Since the critical baryon numbers N crltl 
and N crit2 agree, the number of stable equilibrium solutions is 1 for N < N crit3 , 2 for N crit3 < N < N crit2 = N critl 
and for N crlt2 = N cntl < N. It should be noticed that, for the near critical case —4.9 < (3 < —4.35, the structure of 
the cusp catastrophe does appear although the behavior becomes somewhat more complicated around the maximum 
baryon number for tpo ~ 0. This complicated feature agrees with the fact that the critical mass against zero-mode 
instability is not a monotonic function with respect to /3, which is seen in Table I of p3] . 

Here we comment on the continuous change of the asymptotic value of the scalar field (po . If wc identify i pg w ith 
the cosmological value of the scalar field, the evolution of (fo can be described by the equation of motion (2.4) in 
the Friedmann-Robertson- Walker universe. On the other hand, if wc identify tpo with the value of the scalar field at 
the matching region in the iV-compact-body problem, ipo should evolve due to the change of the density distribution 
around the neutron star. If the time scale of the variation of <pa is sufficiently longer than that of the local gravitational 
phenomena, such as the scalar gravitational wave radiation, the process due to the change of <po can be regarded as 
quasistatic. Through the cosmological evolution of the scalar field ipo, the neutron stars may collapse and radiate a 
scalar gravitational wave. 

We also comment on the continuous change of N, which may be a result of a mass accretion onto the neutron 
star. If the baryon number of the neutron star exceeds the maximum value, the neutron star collapses and scalar 
gravitational waves are radiated and this is a candidate for the source of the scalar gravitational waves [^2| [l4j . In a 
theory like the one of Fig. 7, there is a stable general relativistic neutron star that has the same baryon number and 
ADM energy within numerical accuracy as the maximum-mass non-general-relativistic neutron star has. Then, the 
transition of the non-general-relativistic neutron star to the general relativistic one due to a mass accretion occurs 
without any energy extraction. 

Scalar-tensor theories of gravity naturally arise from the low-energy limit of string theory or other unified theories. 
For the moment, however, it is not clear how the scalar fields should couple to gravity (but see |^5|). Experimental 
tests, such as binary pulsar timing observations, may constrain the way of coupling between the scalar fields and 
gravity. In particular, as for the case in which the single, massless scalar field couples to gravity with the coupling 
function A(cp) = exp[(l/2)/3( j c 2 ], Damour and Esposito-Farese jl6| obtained the constraint on j3 as (3 > —5, using the 
data of three binary pulsars. They showed that the occurrence of spontaneous scalarization makes it very difficult 
for the theory to maintain consistency with the results of binary pulsar timing experiments. The results obtained in 
this paper show that spontaneous scalarization is not an exceptional but robust phenomenon for the neutron star and 
common to a wide range of coupling functions. Gravitational experiments with high- precision and/or in a strong-field 
regime and gravitational wave observations may have the potential to constrain the way of coupling of the gravitational 
scalar fields and thereby we may catch a glimpse of string-scale physics. 
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FIGURE CAPTION 



1. A family of curves of p c — const in the (ipo,N) plane for the j3 = —4 case. The ordinate is nibN in place of 
the baryon number N. In region (A), only one stable equilibrium solution exists, while, in region (B), no stable 
equilibrium solution exists. At a point on the envelope of the family of the curves, the fold catastrophe occurs. 

2. Curves ((po,^) with N constant for the (5 = — 4 case. The number attached to each curve is nibN in the solar 
mass unit. The solid lines denote stable branches and the dotted lines denote unstable branches. 

3. Same as Fig. 1, but for f3 = —6 case. In region (A), only one stable equilibrium solution exists. In region (B), 
two distinct stable equilibrium solutions exist. In region (C), no stable solution exists. Point a is a cusp point. 
At point a, the cusp catastrophe occurs, while the fold catastrophe occurs at a point on the envelopes except 
for a. 

4. Same as Fig. 2, but for the /3 = — 6 case. 

5. (a) (p c ,m) and (b) (p c ,N) curves with ip — for (3 — —6. The solid lines denote stable branches, while 
the dotted lines denote unstable branches. The two distinct bifurcated branches are degenerate because they 
have identical ADM energies and baryon numbers but scalar fields of the opposite sign. The number of stable 
equilibrium solutions changes as 1, 2, as A increases. 

6. Equilibrium space M m in the (ipo, N, w) space around cusp point a for the (3 = — 6 case. This structure of the 
equilibrium space is classified as the cusp catastrophe. 

7. Same as Fig. 5(b), but for (3 = —4.5. The number of stable equilibrium solutions changes as 1, 2, 3, 1, as A 
increases. Two solutions are degenerate on the non-general-relativistic branches. 

8. Same as Fig. 5(b), but for /? = —4.85. The number of stable equilibrium solutions changes as 1, 2, 3, 2, as 
N increases. Two solutions are degenerate on the non-general-relativistic branches. 
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